Reflection, Transmission and Trapping Dynamics of Lattice Solitons at Interfaces 



Y. Kominis and K. Hizanidis 

School of Electrical and Computer Engineering, National Technical University of Athens,Zographou GR-15773, Greece 

Surface soliton formation and lattice soliton dynamics at an interface between two inhomogeneous 
periodic media are studied in terms of an effective particle approach. The global reflection, transmission 
and trapping characteristics are obtained in direct analogy to linear Snell's laws for homogeneous media. 
Interesting dynamics related to soliton power-dependent formation of potential barriers and wells suggest 
spatial filtering functionality of the respective structures. 
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The formation and the properties of spatially localized 
waves in periodic nonlinear structures, known as Lattice 
Solitons (LS), have been the subject of intense research in- 
terest, both theoreticallyand experimentally, in the areas of 
Nonlinear Optics 0,0,13] as well as in Bose Einstein Con- 
densates in optical lattices, [4] while they are also closely 
related to solid state physics. In the context of Nonlinear 
Optics, these studies have led to the emergence of new re- 
search areas where novel phenomena are observed in pho- 
tonic crystals, photonic crystal fibers, coupled waveguide 
arrays and optically induced lattices llLlilj]. Properties of 
wave localization and dynamics in such media with engi- 
neered characteristics suggest potential applications related 
to all-optical signal processing and optical circuits. 

The transverse inhomogeneity of these media results in 
breaking of the translational invariance. Thus, the exis- 
tence of traveling waves is generally not ensured, while 
wave localization can take place in specific positions with 
respect to the underlying structure ||5|,|g]. Depending on the 
complexity of the structure the position and the stability of 
a LS may additionally depend on its power and width (7Q. 
In addition to infinite periodic or quasiperiodic structures, 
LS can be formed at the interface between two semi-infinite 
media with different characteristics. It has been shown that 
such states can be formed at the boundaries between a pe- 
riodic and a homogeneous medium (either linear or non- 
linear) Jsl, 0, [lQ. [Till, a s well as between two dissimilar 
periodic media IllL Il2h. These states, known as Surface 
Solitons (SS), are analogous of the so-called Tamm states 
in semiconductors. 

While the existence of stable and unstable stationary SS 
has been extensively studied in a variety of configurations, 
the investigation of LS and SS dynamics has not been ex- 
plored yet. It is well known that, although the breaking of 
translational invariance does not allow for exact traveling 
wave solutions, LS can travel, without significant radiation 
emission, across an inhomogeneous structure under certain 
conditions for which the effective potential seen by a LS is 
weak. This can be either due to a small depth of the mod- 
ulation of the respective property of the medium or due to 
the fact that the modulation period of the lattice is much 
shorter than the LS width. 

The present work is motivated by the interest on the re- 



flection and transmission characteristics of such mobile LS 
at interfaces between different media and the relation of 
their dynamics with the existence of SS which are trapped 
in the corresponding boundaries. We study LS dynamics 
and SS formation at the interface between two dissimilar 
nonlinear periodic media (the case of an interface between 
a periodic and a homogeneous medium can be treated as a 
special case). A theory is developed, where the formation 
of SS as well as the reflection and transmission character- 
istics of a LS incident, at an angle, at the interface (sepa- 
rating two periodic media), is described analytically by an 
effective particle perturbation theory. The effective poten- 
tial for each soliton determines the location and the stabil- 
ity type of each SS, which is different for waves of differ- 
ent power and width. More importantly, this potential de- 
termines the LS transmission characteristics in an analogy 
to Snell's law of geometric optics. Interesting nonlinear 
dynamics related to the presence of potential barriers and 
wells are demonstrated. The dynamical features are shown 
to depend strongly on the soliton characteristics, so that 
different solitons propagating at the same medium can pos- 
sess qualitatively different dynamical behavior, suggesting 
a power-dependent spatial filtering mechanism. 

The underlying model describing wave propagation in 
an inhomogeneous, Kerr- type nonlinear medium is the per- 
turbed Nonlinear Schrodinger Equation (pNLSE) 

l °'u — en (x)u = (1) 
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where z and x are the normalized propagation distance and 
transverse coordinate, respectively, e is a small dimension- 
less perturbation parameter, and n (x) is the transversely 
inhomogeneous linear refractive index of the medium. The 
latter describes a structure consisting of two semi-infinite 
lattices interfaced at x = and having, in general, different 
amplitudes and periods 

n (x) = A(x) sin [K(x)x + (j)] (2) 

with (j) a constant phase, A(x) = + (AA/2)[1 + 
tanh(ax)] and if (x) = + (AK/2)[1 + tanh(ax)]. 
The amplitude and the wavenumber of the left side lat- 
tice are (A^~\ K^), while for the right side lattice we 
have = (^ (_) + AA,K^ + AK). The 

smoothness of the transition at the interface (at x = 0) is 
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determined by the parameter a. Note that for = or 
= we have the case of a homogeneous nonlinear 
medium interfaced with a nonlinear periodic lattice. 

For small perturbations (small e) the soliton is treated as 
an effective particle 11311 . whose center of mass moves in 
an effective potential which determines the energy "land- 
scape" seen by each soliton. In order to ensure robust soli- 
ton propagation close to the interface, we consider cases 
where the amplitude and wavenumber differences (A^4 and 
AK) are not very large. In the opposite case, a stable soli- 
ton of one lattice emits a significant amount of radiation as 
approaches the interface in order to transform to a stable 
soliton of the other lattice (or a surface mode) resulting in 
dissipative dynamical behavior. 

Following the standard effective particle approach fl3h 
the motion of the center of mass x is given by the simple 
Newtonian equation of motion 

az z ox 

with m — J \u\ 2 dx is the effective soliton mass and the 
effective potential V e f f is given by 

/•+oo 

V eff (x ) = 2 n (x)\u(x - x )\ 2 dx (4) 



where for the evaluation of the effective potential one 
can use the solution of the unperturbed (e = 0) NLSE: 
u(z, x) — ?7sech [rj(x — x )] e lvx / 2 + 2a ? with dx /dz = v 
and da/dz = —v 2 /8 + r] 2 /2. The effective particle ap- 
proach is valid when the initial wave profile (in the form of 
a soliton) does not break up into multiple solitons and/or ra- 
diation, and this is the case we are interested in this work. 
The resulting effective potential depends strongly on the 
relation between three characteristic spatial scales, namely, 
the soliton width (~ r? -1 ), the lattice period (~ (K ± ) _1 ) 
and the smoothness of the properties variation at the in- 
terface (~ a -1 ). The general form of the linear refrac- 
tive index profile n (x) © does not lead to an analytical 
derivation of the respective effective potential from Eq. ©. 
However, the limiting form of V e f f far from the interface 
(at x G — > ±oo) can be obtained analytically and it is iden- 
tical to the effective potential of a single lattice 

lim V eff (x ) = V {±) sin(K (±) x + (j)) (5) 

withF (±) = e^^Trif^J/sinh [ttK^ /(2rj)]. The 
properties of the formation and the dynamics of solitons 
can be directly obtained from the constant Hamiltonian 
(energy) of the equation of motion © H = mv 2 /2 + 
Vef f ( x o) from which the phase space of the respective sys- 
tem can be readily obtained. 

Far from the interface, similarly to the case of a single 
lattice, stationary LS are formed at the extrema of the po- 
tential with LS with x at the minima (maxima) of the po- 
tential corresponding to stable (unstable) modes, for soli- 
tons with any 77 |5|,LZ|]. While the asymptotic form of the po- 
tential is identical for all 77, its depth depends strongly 



on the ratio K^ ±>} /rj. This depth determines the critical 
value of the soliton velocity v c = (2V^ /m) 1 ^ 2 above 
which a stable LS can be de-trapped from the potential 
minimum and move across the lattice. 

On the other hand, close to the interface not only the 
strength but also the form of the potential depends strongly 
on the characteristics of the soliton, namely on 77. A variety 
of qualitatively different forms of the effective potentials 
can occur, in the same structure for different solitons, as 
we show in the following. The number and the position of 
local minima and maxima of the effective potential provide 
important information for both the statics and the dynam- 
ics of solitons in such structures, that is the location and 
stability of stationary SS as well as the transmission and 
reflection of mobile LS. 
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FIG. 1: (Top) Effective potential (red online) and no(x) (dot- 
ted line, out of scale), and (Bottom) phase plane for an effec- 
tive particle corresponding to a soliton with 77 = 0.6 (Left), 
77 = 0.3 (Right). Trapped (blue online), reflected (red online) 
and transmitted (green online) orbits are shown. The underlying 
structure parameters are e = 0.01, = 0.5, — 0.5, 

if(-) = #(+) = 1, a = 1, 6 = 0. 



Firstly, we consider a structure where the interfaced lat- 
tices have the same wavenumber = = 1 and 
different amplitudes A (_) = 0.5,A (+) = 1. In all the 
following examples, we take e = 0.01 in order to ensure 
validity of the effective particle approach, and we consider 
a smooth property variation having a = 1. For the case 
where = 0, the effective potential is shown in Figs. 
l(top) for two solitons of different power (and width) hav- 
ing 77 = 0.6 and 77 = 0.3. The respective phase plane 
diagrams (x , v) can be directly obtained as in Figs. 1 (bot- 
tom), so that both soliton statics, stability and dynamics can 
be determined. This diagram plays the role of a nonlinear 
analogue of Snell's law, describing LS transmission and 
reflection. Stable (unstable) LS or SS correspond to min- 
ima (maxima) of the effective potential. Close to the stable 
stationary LS or SS, trapped states with v 7^ can oscil- 
late in the respective potential well. While stationary LS 
formed quite far from the interface are trivially located in 
the refractive index minima (stable) or maxima (unstable) 
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FIG. 2: Propagation of a LS with 77 = 0.6 (Left), 77 = 0.3 
(Right) with initial position xo = —17n/2 (Left) , xo = 15tv/2 
(Right), corresponding to minimum of V e ff, and initial velocity 
v = -14 x 10~ 2 (Left, Top), v = -16 x 10~ 2 (Left, Bottom), 
v = 7.2 x 10~ 2 (Right, Top), v = 8.4 x 10~ 2 (Right, Bottom). 
The parameters of the structure are as in Fig. 1. 



for every soliton width 77, the situation is drastically differ- 
ent in the area close to the interface. For 77 = 0.6 [Fig. 
l(left)], the pattern of the effective potential is similar to 
the underlying linear refractive index profile. LS reflection 
can occur for appropriate v only for waves incident from 
left to right with initial v larger but lower than v^~\ 
while LS with initial velocities larger than are trans- 
mitted through the interface [Figs. 2(left)]. The situation 
is qualitatively different for solitons having 77 = 0.3 [Fig. 
1 (right)], where an effective potential barrier, associated to 
an unstable SS, can additionally lead to reflection of waves 
incident from the right hand side [Figs. 2(right)]. LS trans- 
mission through the interface requires higher v due to the 
interface induced barrier. The formation of this potential 
barrier cannot be attributed to either of the two lattices but 
to their interfacing, and it does not have a uniform effect 
on all solitons. 
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FIG. 3: Effective potential (red online) and no (x) (dotted line, out 
of scale) for a LS having 77 = 0.3 in a structure with parameters 
as in Fig. 1, except that cj) — tt/2 (Left), cj> — tt (Right). 



The constant phase (j) [related to the amplitude of the 
linear refractive index at the interface (at x = 0)] plays a 
nontrivial role in the form of the effective potential and the 
corresponding soliton statics and dynamics in the structure. 
For (j) — 7r/2 [Fig. 3(left)] and for solitons with 77 = 0.3, 
we have two local minima of the effective potential in both 
sides of the origin corresponding to stable SS. The case of 




FIG. 4: Propagation of SS with r\ — 0.3, initial position xo = 1.2 
[global minimum of V e ff m Fig- 3 (left)] and initial velocity v = 
5.4 x 10 (_2) (Left), v = 6.6 x 10 (_2) (Right). 



(j) — tt and a soliton with 77 = 0.3 [Fig. 3(left)], results to 
an effective potential of opposite sign with the one shown 
in Fig. l(top, left). The global minimum of the potential 
close to the interface corresponds to a stable SS. Solitons 
with initial velocities, corresponding to a total energy level 
inside the potential well are trapped and oscillate around 
the position of the potential minimum [14]. The corre- 
sponding de-trapping velocities can be directly obtained 
from the effective potential. This is an important advan- 
tage of the effective particle approach because, in addition 
to determining the stability of a stationary SS (as linear sta- 
bility analysis could also do), it provides a measure of the 
degree of stability of the SS in terms of the depth of the 
effective potential and the de-trapping velocity. The latter 
is directly related to the launching angle of the soliton at 
z = 0, and gives the tolerance on the alignment precision 
for soliton trapping. In Fig. 4(left), we show a trapped 
SS with initial energy less than the local maximum located 
on the left of origin, oscillating around the position of the 
potential minimum. For velocities corresponding to initial 
energies larger than the local potential maximum on the left 
of the origin [Fig. 3(left)], the amplitude of the soliton po- 
sition oscillation is larger [Fig. 4(left)], as predicted from 
the effective potential. 
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FIG. 5: Effective potential (red online) and no (x) (dotted line, out 
of scale) for a wide LS with 77 = 0.1. The underlying structure 
parameters are as in previous Figs, except that = (Left), 
cj) = tt/2 (Right). 

From the previous cases it is shown that the formation 
of nontrivial effective potential patterns including barriers 
and wells, in the vicinity of the interface, is mostly related 
to solitons with large (in comparison to the lattice period) 
widths (77). It is remarkable that for wide LS, although 
the amplitude of the effective potential is very small suf- 
ficiently far from the interface [as obtained by Eq. ©], it 
has large values close to the interface. This means that al- 
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though a wide LS does not "feel" the lattice, it does "feel" 
the interface. Such cases, for a LS with r\ = 0.1, are 
shown in Figs. 5 for = 0(left) and <j) — tt/2. Ac- 
cording to the respective effective potentials, a LS which 
is untrapped from the lattice far from the interface and can 
move, it can be reflected from the interface even if it has a 
large v [Fig. 5 (left)]. Also, stable and trapped SS can be 
formed in the vicinity of the interface in a position depend- 
ing on the phase constant (j) [Fig. 5(right)]. 
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FIG. 6: Effective potential (red online) and no (x) (dotted line, out 
of scale) for LS having 77 = 0.9 (Left), 77 = 0.4 (Right) in a struc- 
ture consisting of two interfaced lattices with different periods. 
Parameters: e = 0.01, = 0.5, = 0.5, if (_) = 0.5, 

=1, a = 1,0 = 0. 




pie in Figs. 7. This difference is manifested as a displace- 
ment of the reflected wave from its geometric optics path. 
For the case of linear wavepackets, this effect is referred as 
the Goos-Hancen shift. For the nonlinear case of LS soliton 
reflection at an interface and within the effective particle 
approximation the natural analogue of the Goos-Hancen 
shift 5z is defined as the difference in the asymptotic value 
of the z location of the effective particle as x — > ±00 
and the location of a free particle bouncing at the inter- 
face x = B13I1 . The respective analytical expression is 
5z = J^ r 00 [v~ 1 (x) — VQ^dx + x r /v , where v 2 (x) = 
Vq — (2/ m) [V e f f (x) — V e f f (±00)] and x r is the reflection 
position given by v% = (2/m)[V eff (x r ) - V eff (±oo)]. 

In conclusion, by utilizing an effective particle approach, 
shown in good agreement with propagation simulation re- 
sults, we have studied SS formation and LS dynamics at 
the interface between two dissimilar inhomogeneous me- 
dia. The global trapping, reflection and transmission char- 
acteristics have been described in a direct analogy to Snell' 
laws for linear homogeneous media, while an analytic ex- 
pression for the nonlinear analogue of the linear Goos- 
Hanchen shift has been obtained. It has been shown that 
different solitons can have qualitatively different dynam- 
ics in the same structure due to power (or width) depen- 
dent formation of effective potential barriers and wells. 
These effects could conceptually be considered as a power- 
dependent spatial filtering functionality. 



FIG. 7: Reflection of a LS with 77 = 0.9 (Left), 77 = 0.4 (Right) 
incident at the interface from an initial position xq = —IStv/2 
(Left), x = 15tt/2 (Right) with v = -19.6 x 10~ 2 (Left), v = 
12 x 10~ 2 (Right). 

Finally, an interesting case occurs when not only the 
amplitudes of the two interfaced lattices are different, but 
also the wavenumbers (K^). In Figs. 6, we show 
the corresponding effective potentials for solitons having 
77 = 0.9 (left) and 77 = 0.4 (right) in a configuration where 
A<-> = 0.5, #' = 1, = 0.5 and = 1. It is 
interesting that the same structure results in a qualitatively 
different "landscape" for different solitons. Therefore, for 
solitons with 77 = 0.9, < and reflection can 

take place only for solitons incident from the left, while 
the opposite holds for solitons with 77 = 0.4 (Figs. 7). For 
solitons with 77 ~ 0.6 no reflection can take place. There- 
fore, a power-dependent directional effect is introduced. 

In the previous cases we have investigated, in terms of 
the effective potential (and the corresponding phase plane), 
whether a soliton with a given initial position and velocity 
can be trapped, reflected or transmitted in the interface be- 
tween two dissimilar lattices. In the case of reflection, an 
additional information of interest is related to the position 
where a LS is actually reflected, which can be quite differ- 
ent under certain parameter selections as shown for exam- 
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